Demonstration of unconditional one-way quantum computations for continuous 

variables 
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Quantum computing promises to exploit the laws of quantum mechanics for processing informa- 
tion in ways fundamentally different from today's classical computers, leading to unprecedented 
efficiency [l], Q|. One-way quantum computation, sometimes referred to as the cluster model of 
quantum computation, is a very promising approach to fulfil the capabilities of quantum informa- 
tion processing. The cluster model is realizable through measurements on a highly entangled cluster 
state with no need for controlled unitary evolutions [H, [3- Here we demonstrate unconditional 
one-way quantum computation experiments for continuous variables using a linear cluster state of 
four entangled optical modes. We implement an important set of quantum operations, linear trans- 
formations, in the optical phase space through one-way computation. Though not sufficient, these 
are necessary for universal quantum computation over continuous variables, and in our scheme, in 
principle, any such linear transformation can be unconditionally and deterministically applied to 
arbitrary single-mode quantum states. 



The cluster model of quantum computation (QC) is 
a recently proposed alternative to the conventional cir- 
cuit model [3j-|8|. In this model, unitary operations are 
achieved indirectly through measurements on a highly 
entangled quantum state - the cluster state. Cluster 
computation is achieved through the following steps: (1) 
preparation of an entangled cluster state and an input 
state for processing; (2) entangling operation on these 
two states; (3) measurements on most subsystems of the 
cluster state and feed-forward of their outcomes; (4) oc- 
currence and read-out of the output in the remaining 
unmeasured subsystems of the cluster. Universality, i.e., 
realization of arbitrary unitary operations is achieved by 
adjusting the measurement bases, sometimes also depen- 
dent on the results of earlier measurements [f| . 

Several experiments of one-way quantum computation 
have been reported for discrete-variable (qubit) systems 
using single photons These demonstrations of 

one-way quantum computation work in a probabilistic 
way, since the resource cluster is generated only when 
the photons that compose the cluster are produced and 
detected. Another typical feature of the single-photon- 
based cluster computation experiments is that the usual 
input states, |+), are prepared as part of the initial clus- 
ter states. These properties would pose severe limitations 
when unitary gates are to be deterministically applied on- 
line to an unknown input state which is prepared inde- 
pendently of the cluster state, for instance, as the output 
of a preceding computation. 



In contrast, we report in this paper on unconditional 
one-way quantum computation experiments conducted 
on independently prepared input states. These inputs, 
as well as the entangled cluster state, are continuous- 
variable states. The price to pay for this is a set of 
stronger requirements on universality. Not only do we 
need at least one nonlinear element to achieve completely 
universal QC over continuous variables 0, we also 
have to cover all linear transformations, which, for a sin- 
gle optical mode, consist of arbitrary displacement, ro- 
tation, and squeezing operations in phase space. Our 
scheme represents the ultimate module for arbitrary lin- 
ear transformations of arbitrary one-mode quantum op- 
tical states. It can be directly incorporated into a full, 
universal cluster-based QC together with a nonlinear el- 
ement such as measurements based on photon counting 
[3| (for a discussion on the fidelity when concatenating 
our module using finitely squeezed cluster states, and 
on its scalability into a full, measurement-based QC, see 
supplementary information and Refs. 15-17j|). 

We use a continuous-variable four-mode linear clus- 
ter state as a resource Q. An approximate version 
of this cluster state can be obtained deterministically 
by combining four squeezed vacuum states on an 80%- 
transmittance beam splitter and two half beam splitters 
(HBSs) [Hi^. 

Recently, it was shown that the complete set of one- 
mode linear unitary Bogoliubov (LUBO) transforma- 
tions, corresponding to Hamiltonians quadratic in x and 
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FIG. 1: (A) Abstract illustration and (C) experimental setup of one-mode LUBO transformations using a four-mode linear 
cluster state. There is a 1-to-l correspondence between (A) and (C). Squeezed vacuum states are generated by subthreshold 
optical parametric oscillators containing periodically poled KTiOPCU crystals as nonlinear media. (B) phase space represen- 
tations of quantum states in each step of the Fourier transformation (B-i) and the lOdB rr-squeezing operation (B-ii), starting 
with a vacuum state input (blue) and an ^-coherent state input (red). Tele.: teleportation, Op: operation. 



p, can be implemented using a four-mode linear cluster 
state as a resource (2lj . The measurements required to 
achieve these operations are efficient homodyne detec- 
tions with quadrature angles 9i, which are easily control- 
lable by adjusting the local oscillator phases in the homo- 
dyne detectors. The total procedure then consists of the 



tclcportation-based [22H24| coupling M te i e (9i n ,9i), fol- 



lowed by two elementary, measurement-based, one-mode 
operations M(9i 
mation): 



14J . |25| . |26| (see supplementary infor- 



IVw) = M(9 3 )M(9 2 )M tele (9 m ,9 1 )\ij in ) 



(1) 



Each step can be decomposed into three inner steps, 
a ^-rotation, squeezing, and a ^-rotation in phase 
space: R(<p)S(r)R(</)) with R(9) = e w ^ 2+ ^ and 
S(r) = e^P+^l We have M tele {e in A) = 

R(-9 + /2)S(r)R(-9 + /2) with r = logtan(0_/2) and 
6± = 9 m ± 6 X , while M{6i) = R^S^R^i) with 
Ti = log V k i+ tt^l^ ( f >i = I — tan -1 V fc i+ 4 k i _ an j fc. _ 
1/ tan 9i. 

In our experiment, we demonstrate four types of LUBO 
transformations: the Fourier transformation F = R(%/2) 
(90° rotation) ; and three different x-squcczing operations 
S{r) with r = In 10*, a = 3,6,10[dB]. FIG. UK and 
FIG. [Tp show the abstract illustration and the experi- 
mental setup, respectively. We employ the experimental 
techniques described in Refs. [IH and for the gener- 
ation of the cluster state and the feed-forward process, 
respectively. 

The Fourier transformation is achieved by choosing for 
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FIG. 2: Fourier transformation operation; (A) Measurement 
results of the input state. Trace (A-i) shows the shot noise 
level (SNL) and (A-ii) shows the phase scan of the input 
state. (B) Measurement results of the output state. Trace 
(B-i) shows the SNL, (B-ii) shows the phase scan of the out- 
put state, and (B-iii) shows the measurement result of the x 
quadrature with a vacuum input. The measurement quadra- 
ture angle is determined through the relative phase between 
the signal beam and the local oscillator beam. The measure- 
ment frequency is 1 MHz and the resolution and video band- 
widths are 30 kHz and 300 Hz, respectively. Traces (A-i), 
(B-i), and (B-iii) are averaged 20 times. 



step (3) measurement quadrature angles (9i n , 61, 62, 03) 
as (90°, 0°, 90°, 90°), see supplementary information. 

The measurement results for the Fourier transforma- 
tion of a coherent state input are shown in FIG. [5J As 
clearly shown in FIG. HJA., the input is a coherent state 
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Time (s) Tim e (s) Squeezing (Operation) (dB) 

FIG. 3: Squeezing operations; (A, B) lOdB ^-squeezing operation with an a;-coherent input(A) and a vacuum input(B). (i) 
shot noise level; (ii) phase scan of the output state; (iii) measurement of x; and (iv) measurement of p. The measurement 
settings are the same as in FIG. [2] Traces (i), (iii), and (iv) are averaged 20 times. (C) experimental results (dots) and 
theoretical calculation (solid curves) of 3dB, 6dB, and lOdB x-squeezing operations. Black and blue traces correspond to a p 
measurement with p-coherent input, and x measurement with :r-coherent input, respectively; green and red traces correspond 
to a p measurement with vacuum input, and x measurement with vacuum input, respectively. Each data point has an error of 
about ±0.2dB. 



with amplitude 17.7±0.2dB. The output state is shown in 
FIG. dp. The peak level of trace FIG.[2£(ii) is 17.5±0.2 
dB higher than the shot noise level (SNL), which is the 
same level as the input within the error bar. We acquire 
the peak of the input by measuring x, while we obtain 
the peak of the output by measuring p, corresponding to 
a 90° rotation in phase space. These measurement re- 
sults confirm that the Fourier transformation is applied 
to the input coherent state. 

The quality of the operation can be quantified by 
using the fidelity, defined as F — (\P Ueal\Pout\^ 'ideal) ■ 
In the specific case of our experiment, the fidelity 
for a coherent input state as given above is F = 
2Ml+4< tt )(l+4< rf ), 



where a x out 



and <7„„, t are the 



variances of the position and momentum operators in the 

Wc obtain a* 



output state, respectively [31 



2.9±0.2 



dB(FIG. and a p out = 2.8±0.2 dB (not shown) 

above the SNL with a vacuum input, corresponding to 
a fidelity of F=0.68 ± 0.02. This is in good agreement 
with the theoretical result F = 0.71, where an average 
squeezing level of — 5.5dB is taken into account. 

Another fundamental element of the LUBO trans- 
formations is squeezing. A sequence of teleportation 
coupling M te i e (0i n ,9i) followed by elementary one-mode 
one-way operations M{9{) is required in order to extract 
squeezing without rotations (see FIG. [Tj3(ii)). 

Wc implemented three different squeezing operations 
with three different sets of quadrature measurement an- 

gics {e in ,6 1 ,e 2 ,e 3 y. 



(-42.5°, 62.4°, 63.5°, 76.0°) 
(-41.4°, 72.2°, 41.9°, 74.4°) 
and (-47.7°, 79.2°, 25.9°, 78.4°) 



(2) 



resulting in 3dB, 6dB, and lOdB x-squeezing opera- 
tions, respectively (see supplementary information). In 



all these squeezing gates, the inputs are chosen to be co- 
herent states with a nonzero amplitude in x (x-coherent) 
or in p (p-coherent), and these amplitudes are 14.7dB± 
0.2dB. 

FIG. [3]4. shows the measurement results of the lOdB 
x-squeczing operation on the x-coherent state. In this 
figure, the extra dotted lines are plotted for compari- 
son, in order to show the levels of the input state: x in 
blue (14.7dB) and p in black (SNL). Wc obtain signal 
levels of 5.1±0.2dB and 11.5±0.2dB above the SNL for 
the measurement of the x and p quadratures of the out- 
put, respectively. The level of the x quadrature of the 
output (FIG. EK(iii)) is about lOdB lower than that of 
the input (the blue dotted line in FIG. OV), while the 
variance of the p quadrature of the output (FIG.[3J^.(iv)) 
increases by about lOdB compared to that of the input 
(the black dotted line in FIG. [3]A.). These observations 
arc consistent with a lOdB x-squeezing operation. Note 
that the x and p quadratures of the output have addi- 
tional noises. These are caused by the finite squeezing of 
the cluster state and would vanish in the limit of infinite 
cluster squeezing. 

In order to show the nonclassical nature of the out- 
put state, we also use a vacuum state as the input 
(FIG. |3)3). The measured variance of the x quadrature 
is — 0.5±0.2dB, which is below the SNL, thus confirming 
nonclassicality. 

Finally, we demonstrate the controllability of the one- 
way quantum computations. Both theoretical curves 
(with ~5.5dB resources) and measured results for the 
three levels (3dB, 6dB, and lOdB) of x-squeezing are plot- 
ted in FIG.|3P. Three kinds of input states are used here: 
a vacuum state; an x-cohercnt state; and a p-coherent 
state. As can be seen in FIG. |3P, the measurement 
results agree well with the theoretical curves, and all 
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the operations are indeed controlled by the measurement 
bases for the four homodyne detections. 

In summary, we have experimentally demonstrated 
one-way quantum computations with continuous vari- 
ables. All operations were perfectly controllable through 
appropriate choice of measurement bases for the homo- 
dyne detections. In our scheme, arbitrary linear one- 
mode transformations can be applied to arbitrary input 
states coming independently from the outside. An exten- 
sion to multi-mode transformations, though not demon- 
strated here, is also possible by similar means [21|. The 
accuracy of our one-way quantum computations only de- 
pends on the squeezing levels used to create the resource 
cluster state. Although in our experiment squeezing lev- 
els were sufficient to verify the nonclassical nature of 
the output states, even higher levels of squeezing, as re- 
ported recently [32|, 



may lead to increased accuracies 
and one-way quantum computations of potentially larger 
size in the near future. In order to achieve quantum 
operations other than linear unitary mode transforma- 
tions, nonlinear measurements besides homodyne detec- 
tions would be required. However, the demonstration of 
the experimental capability of implementing an arbitrary 
linear single-mode transformation through continuous- 
variable cluster states, as presented here, represents a 
crucial step toward universal one-way quantum compu- 
tation. 



Appendix B: Continuous- Variable Cluster 
Computations 

In this experiment, we use a continuous-variable four- 
mode linear cluster state: 

dadbe 2iab \ Pl = a)\x2 = a)\x 3 = b)\p 4 = b), (5) 

as a resource, where \x = c) and \p = c) = F\x = c) 
(with the Fourier transformation F = e^^ x + v >) are 
cigenstates of the canonical conjugate position and mo- 
mentum operators, respectively x andp, with eigenvalues 
c € R (fi = 1/2); the subscripts label the corresponding 
modes. Here, \x) is the computational basis for our CV 
system. 

As is mentioned in the main text, the output state 
becomes 



hkut) = M{e z )M(e 2 )M tele {e in ,e 1 )\^ v 



(6) 



Note that M te i e (0,- n , 6\ ) cannot be decomposed into 
Mi(6i)Mi n (9i n ), because the measurements on modes 
in and 1 are nonlocal measurements. The operations 
Mtei e {@ini@i) an d M{0i) are each elements of the one- 
mode LUBO transformations. 



In the following sections, we show explicit derivation 
of these operations and measurement quadratic angles. 



Appendix A: Discrete- Variable Cluster 
Computations 



Appendix C: Quantum Computation using a 
Four-mode Linear Cluster State 



In the experiments reported in Rcfs. quantum 
computations are demonstrated by showing arbitrary ro- 
tations on a qubit using a discrete-variable four-qubit 
linear cluster state: 



E 

{a,6} = {0,l} 



(-ina)|a)|fe)|&), 



(3) 



where |0) and |1) are the computational basis states, 
while |0) = |+) = (|0) + |1))/V2 and |1) = |-) = 
(|0) — \l))/^/2 can be obtained with |i) = H\i) and the 
Hadamard gate H. Using this cluster state as a resource 
for cluster computation (note that, in the recent singlc- 
photon-based works, the usual input states, |+), are pre- 
pared as part of the initial cluster states), a sequence of 
four operations can be applied onto an input state \ip)i n , 



i 



(4) 



e -i4>iZ/2 are 



where = e -^* x / 2 and Z^ 

rotations about the X and Z axes on the qubit's Bloch 
sphere with the usual Pauli operators X and Z, respec- 
tively. 



In the Heisenberg picture, a perfect four-mode linear 
cluster state has zero-eigenvalue correlations: 



PS 

0§ 



P? 



ffi-xg = S 



5x^0 

= S 2 ->0 

6 3 ->-0 

4 ->0, 



(7) 



in the limit of infinite squeezing. The Xj and pj are po- 
sition and momentum operators for an optical mode j 



with an annihilation operator a,- 



tpj. 



In the ex- 



periments, squeezing levels are limited, thus Si have non- 
zero variances. An approximate four-mode linear cluster 
state can be generated by combining four squeezed vac- 
uum states on a 80%-transmittance beam splitter and 



two half beam splitters (HBSs) 
extra noise terms 



[lj, I2Q|, leading to the 



' Si 


= V2e- 


r p S i 






5.2 




- r p§ + 




- r pf 
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1 

= 7f 


- r pf- 




- T P S 2 




= V2e- 









(8) 
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where = e r £j + ie~ r pj shows a squeezing resource 
for the cluster state. We assume that each resource has 
a same squeezing level r. 

Also in our scheme, an unknown input can be coupled 
with the four-mode linear cluster resource using a half 
beam splitter utilizing the process of quantum teleporta- 
tion. The input coupling is expressed as follows: 



■I 1 




(9) 



ipi n shows an annihilation operator 



where a% n = £i n 
for an input mode. 

The modes in, 1-3 are measured simultaneously by us- 
ing homodyne detections with measurement quadrature 
angles 9j, and all feed- forward processes are postponed 
until the end of the cluster computation 21 , [25| . The 
measurement variables are 



x l 

iM 



( £?„ COS 



m 

&i cos 9 1 

£2 COS 62 



■p?sm6 in \ 
Pi sin 61 
p% sin 9 2 



(10) 



\ £3 cos 9 3 + pC sin 6 3 / 



By using the equations above, the x and p quadratures 
of mode 4 can be expressed as 



1 

02 



--M(k 3 )M(k 2 )M tele (6 



M(k 3 )M(k 2 )M M (6^ 



M{k 3 ) 



xf /sin ; 







M{k 3 )M{k 2 



M(k 3 ) 



where 



M M (0+,0-) = 



-1 
1 



64 



' cos 9- -f-cos 9+ 
sin 6- 
— sin 9+ 



Pi 

x%/sm6 in 
£f /sin 0i 

if/sin0 3 ' 
Si 



(11) 



sin y_ 

cos 9+ — cos ( 



\ sin 9- sin 9- 

sin 9- +sin sin 9- — sin 9+ ' 
V^sintf- V2 sin 9- 

cos 9— —cos 0+ cos 9- —cos #_|_ 



sin 9- 

l 



>± — V'in 



±0i. (12) 



The first term on the right-hand side of Eq. (TTT|) is the 
main operation controlled by measurement quadrature 
angles. M(ki) and M te;e (0 + ,0„) in Eq. dTTJ with Eq. 
(fl"2]) correspond to M(0 4 ) and M te i e (0 m ,0i) in Eq. ©, 
respectively. The second to fourth terms correspond to 
back-actions of the measurements; the quadrature opera- 
tors of mode 4 are shifted depending on the measurement 



results Xj , and these terms should be eliminated by a 
succeeding feed-forward process. The remaining terms 
show additional components caused by imperfec- 
tion of squeezing resources, which lead to errors in cluster 
computations. 

In our case with Eq. (|5J|, the errors are 



-= - V2k 2 h 



I ^ -r ~S 



5 P = y/2k 2 e~ r pf - \l h~ r p§ + -U^pf, 



I ~j= ~ VZ/t 2 k 3 je -pi - y -e 

thus additional variances are 
((AS X ) 2 ) =±e- 2 rU±-V2k 2 k 3 ^ +l + 3kl 
((AS p f) =ie- 2 '"(3 + 2fcl). 



(13) 



(14) 

Note that Eq. (|13p and Eq. (TT4"]) consist only of squeezing 
components derived from the squeezing resources because 
all antisqueezing components are eliminated in Si (see Eq. 
(|8|)). In this case all antisqueezing components derived 
from the squeezing resources in the output vanish via 
the feed-forward process. This leads to a considerable 
reduction of extra noise terms in the output. 



Appendix D: Derivation of Measurement Angles for 
the Fourier Transformation 



In order to realize the Fourier transformation 



Xout\ _ f0 -l\ ( Xin 
Pout) U \p m l ' 



the following equation should be satisfied: 

M(k 3 )M(k 2 )M tele (6 + ,6-)= (J ~ 2 



Therefore, we get 



kin = 0, k 2 = —,h = 0, 

h 



(15) 



(16) 



(17) 



where ki is a free parameter which can be chosen such 
that the error in the output is minimized. Here, as can 
be seen from Eq. (Ti"4"]) . the additional variances are 

((A5 x f) = \e~ 2r , ((A4) 2 ) = V 2r (3 + 2kj) , (18) 

and ((A(5 p ) 2 ) becomes minimal when k 2 = 0. Thus, 
{kin-,k\,k 2 ,k 3 ) = (0,oo,0,0) is the optimal set, which 
corresponds to (6 in ,9i,6 2 ,6 3 ) = (90°, 0°, 90°, 90°). 
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The relation between the input and output is given by 



Xout\ _ (0 -1 
Pout J V 1 



+ 



/ 1 -r-S 

5 



1 



-e p 3+ ^=e p A 



(19) 



and thus, the Fourier transformation is indeed achieved 
and the additional variances are 

<(A4) 2 > = \e~ 2r , ((A5 p ) 2 ) = \e~ 2 \ (20) 



Appendix E: Derivation of Measurement Angles for 
the a-squeezing Operations 

Next, we move on to adB ^-squeezing operation: 



Xout\ _ AO 2° \ f x in 
Pout J \ 10& J \p in 



(21) 



In order to achieve this operation, kj should be selected 



fci 



A.'; 



1 + 2 • lO^fc,. 



l + 10*fc in 10-*(l + 10*fci„) 
k 2 = r , As = r • (22) 



fcj„ can be chosen such that the quadrature x of the out- 
put has a minimum additional variance. From Eq. (|14p . 
the optimum fci n is independent of the levels of squeezing 
resources r, and it is determined only by the operation 
level a. Straightforward algebra shows that minimum 
variances occur when fc,„ is 



kin = — 10 20 



2 (-2 • 1055 ± V-3 + 2- 10* 

= : —a (10TO > f) 



(io* < I) 



3 + 2 • IOto 



(23) 



Both ki n+ and km- give us the same variance and fc_ 
is selected for our experiment. Therefore, the angles 

(din, &i, 2 , #3) : 



(-42.5°, 62.4°, 63.5°, 76.0°), 
(-41.4°, 72.2°, 41.9°, 74.4°), 
and (-47.7°, 79.2°, 25.9°, 78.4°), 



(24) 



should be selected for 3dB, 6dB, and lOdB x-squeezing 
operations, respectively. 



Appendix F: Remarks on Fidelity and Scalability 

Fidelity- when our module for arbitrary linear one- 
mode transformations is concatenated as needed for 



larger quantum computations, the noise originating from 
the use of realistic, finitely squeezed cluster states will 
lead to an accumulation of errors and a decreasing fidelity 
for the output state. This effect is inevitable and will oc- 
cur even when all the remaining operations including the 
homodyne measurements are performed with 100% effi- 
ciency. 

The quality of a cluster-based operation starting with 
an initial pure state can be quantified using the fidelity 
defined as 



F = (^ideal\Pout\^ 



ideal) 



(25) 



where \^ ideal) and p ou t are the ideal pure output state 
and the density matrix of the experimental output state, 
respectively. We shall consider the measurement-based 
application of Fourier transformations starting with an 
arbitrary pure Gaussian state, and we obtain 



F 



v/(l + 4<J(l+4<J' 



(26) 



where <J% ut and o v out are the variances of the position 
and momentum operators in the output state, respec- 
tively. Since the excess variances for a Fourier transfor- 
mation, realized through four elementary measurement- 
based steps, are given by Eq. (|2"0"j) . the fidelity becomes 



F 



1 



3 -2r 
2 C 



(27) 



We can now easily extend this discussion to the fi- 
delity for an n-step telcportation (or, more specifically, 
a Fourier transformation, assuming n is even). In this 
case, we have an input coupling through telcportation 
(two steps) followed by an (n — 2)-step one-mode one- 
way gate with the fidelity, 



1 



1+ *±l e -2r' 



F = 



fe+l c -2r 



fc+2 -2r 
2 e 



n = 2k, 
-., n = 2A;+1. 



(28) 

Note that a cluster-based n-step quantum teleportation 
roughly corresponds to an ^-step sequential quantum 
telcportation, because the fidelity for an m-step sequen- 
tial quantum teleportation is F — 1 _i_ TT ^ e -a T - ■ The ex- 
cess noise of one-way QC through a linear chain is now 
roughly given by ne~ 2r . Therefore, for a larger computa- 
tion with an increasing n value, but an unchanged output 
fidelity, the required level of resource squeezing (e +2r ) is 
roughly proportional to n. In other words, using the vari- 
as a figure of merit (the so-called accuracy of 



ance e 



+2r 



the cluster state [14j), there is a linear (and not an ex- 
ponential) dependence of the required accuracy on the 
length of the computation. Further, any desired accu- 
racy can be achieved for a cluster state of arbitrary size 
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with the same squeezing levels, provided the connectivity 
(the maximum number of nearest neighbors for any mode 
of the cluster) is constant like in the present example of 
a linear chain fl~H ]. Similarly, the entanglement between 
any mode and the remaining modes of the cluster state 
for a given squeezing level only depends on the number of 
links for that single mode and is independent of the size 
of the cluster [171] . However, note that the noise accu- 
mulated in a computation and hence the accuracy of the 
computation, for a given accuracy of the cluster resource, 
does depend on the length of the computation and hence 
on the size of the cluster, as described above. So there 
is a distinction between the accuracy of the cluster and 
the accuracy of the computation; the former can be in- 
dependent of the size of the cluster, whereas the latter, 
of course, is not. 

Scalability- we may consider the result of cascaded 
quantum teleportations through a linear cluster chain 
with a certain finite error specified by a lower bound Fq 
on the fidelity. In this case, we obtain for the number of 
possible steps, 



4e 



+ 2c 



1 



1 



1. 



(29) 



omitting the parity of n here for simplicity. Now consid- 
ering the "classical limit" of teleportation, Fq = ^ [3l , as 
a benchmark and a squeezing level of about — 5.5dB like 
in our experiment, n = 13 steps of elementary teleporta- 
tions would be possible. A recently reported squeezing 
value of — 12.7dB [35| would enable us to cascade the 
teleportations up to 73 times. This alone shows already 
that in a weak sense (see below), our scheme is scalable 
and can be extended to a higher number of quantum op- 
erations (than just the present four operations) with the 
same experimental parameters and squeezing resources. 

Finally, we shall briefly comment on the full scalability 
of the present experimental scheme, and measurement- 
based (MB) QC over continuous variables in general. An 
undeniable fact is that there is no proof for scalability 
of continuous- variable QC in the presence of errors in a 
strict sense. Strict here means that an analogous result to 
the so-called threshold theorems in the discrete-variable 
regime is still lacking for continuous variables. While it 
has been shown that arbitrarily long qubit computations 
(in a circuit model) can be achieved to any desired ac- 
curacy provided the elementary components of the com- 
putation are less faulty than a certain threshold value 
|l(, such a result does not exist for continuous variables. 
Mapping the circuit-based thresholds to measurement- 
based thresholds is also possible in the discrete-variable 
regime for certain abstract error models [H, |3tJ • How- 
ever, these error models do not directly apply to those 
realistic, experimental situations in which cluster states 
are prepared highly probabilistically through parametric 
down conversion (PDC) by means of linear optical ele- 
ments HHH. 



In the continuous-variable case, even if there was a 
circuit-based threshold theorem, the transition from the 
circuit to the MB model appears to be fundamentally dif- 
ferent compared to the discrete- variable case. The reason 
is that the continuous-variable Gaussian cluster states 
are intrinsically noisy due to their finite squeezings, al- 
ways resulting in squeezing-induced errors in a MBQC 
(see above). In a circuit computation, such errors would 
not occur, except for a possible encoding step into, for 
instance, an approximate position eigenstate. 

There arc now a couple of recent investigations into the 
effects of finite squeezing on the scalability of MBQC with 
Gaussian cluster states. In Ref. [Hj], it is argued that the 
accumulation of squeezing-induced errors prevents scala- 
bility in the strict sense of arbitrarily long computations, 
as long as no extra tools such as quantum error correc- 
tion codes are incorporated from the beginning. More 
precisely, when nonclassical correlations in form of an 
entangled state are to be transmitted through a linear 
continuous-variable cluster chain, the entanglement will 
decay exponentially with the length of the chain, simi- 
lar to the effect of cascaded entanglement swapping with 
two- mode squeezed states [HI . The results of Ref. 15 1 , 
however, are even more general than the simplest case of 
homodync-based swapping along a linear chain, including 
as well non-Gaussian measurements. Contrary to these 
negative observations concerning scalability, it was also 
shown recently that an important necessary requirement 
for a cluster state to be an efficient universal resource [39[ 
can be indeed satisfied by the Gaussian, finitely squeezed 
cluster states [HI]. Similar to this more optimistic view, 
one should also note that, even though the entanglement 
in a linear Gaussian cluster chain does exponentially de- 
cay, there is no need to use an initial resource squeez- 
ing r that grows exponentially with the total number of 
measurement-based computation steps. In fact, the ex- 
ponential decay rate itself depends on r, establishing a 
quantitative link between the necessary initial squeezing 
resources and the maximal number of operations n rnax 
that still result in an effective output squeezing r n ut al 
an arbitrary constant bound c, r out > c, namely [40( 



ou t above 



- e +2r lntanht 



(30) 



Therefore, in particular, we have 2r ~ \x\{n max ) + const., 
corresponding to a logarithmic increase of the initial 
squeezing with the maximal number of operations that 
still satisfies a certain accuracy threshold. Note that this 
statement can be equivalcntly made in terms of an entan- 
glement measure such as the so-called logarithmic nega- 
tivity which for a pure two-mode squeezed state is pro- 
portional to the squeezing parameter. In this case, we 
obtain for the initial resource entanglement E the scal- 
ing property E ~ ln{n max ) + const, to guarantee that 
the output entanglement satisfies E out > c for some ac- 
curacy bound c. The bottom line of our discussion here 
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is that the required input squeezing and entanglement, 
to make sure that the output entanglement along an ar- 
bitrarily long cluster chain remains at least as large as 
some fixed bound, scales logarithmically with the length 
of the chain. 

Finally, compared to the existing theoretical results 
on discrete-variable fault-tolerant QC and the published 
experiments on single-photon-based qubit MBQC, even 
the results of Ref. [15J do not rule out the possibility of a 
fault-tolerant version of MBQC over continuous variables 
(see also the final discussion in Ref. [Hj]). However, in 
order to deal with finite-squeezing errors from the start, 
there is no known error correction scheme shown to be 
capable of suppressing such errors (at least for a logi- 
cal continuous- variable state); it is only known that such 
an encoded scheme must be based upon some nonlin- 
ear, non-Gaussian element j4lj . similar to the nonlinear 
measurement which would be needed to achieve universal 
operations through the cluster state. 
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